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VHE Gamma-rays from Galactic X-ray Binary Systems 
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The detection of TeV gamma-rays from LS 5039 and the binary pulsar PSR B1259-63 by HESS, and from LS I +61 303 
and the stellar-mass black hole Cygnus X-l by MAGIC, provides a clear evidence of very efficient acceleration of parti- 
cles to multi-TeV energies in X-ray binaries. These observations demonstrate the richness of non-thermal phenomena 
in compact galactic objects containing relativistic outflows or winds produced near black holes and neutron stars. I 
review here some of the main observational results on very high energy (VHE) 7-ray emission from X-ray binaries, as 
well as some of the proposed scenarios to explain the production of VHE 7-rays. I put special emphasis on the flare 
TeV emission, suggesting that the flaring activity might be a common phenomena in X-ray binaries. 
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1. Introduction 

An X-ray binary (XB) is a binary system contain- 
ing a compact object, neutron star (NS) or a stellar- 
mass black hole (BH), emitting in X-rays. In most 
cases, the X-ray emission is a consequence of the ac- 
cretion of matter from the companion star. Up to 
now, 299 accreting X-ray binaries have been cata- 
logued in the Galaxy, being 114 High Mass X-ray 
Binaries (HMXBs) [1] and 185 Low Mass X-ray Bi- 
naries (LMXBs) [2]. In the HMXBs, the optical com- 
panion has spectral type O or B and the mass trans- 
fer occurs via decretion disc (Be stars) or via strong 
wind or Roche- lobe overflow. In the LMXBs, the op- 
tical companion has spectral type later than B and 
the mass transfer is via Roche-lobe overflow. Among 
the 299 known XBs, 65 of them (9 HMXB and 56 
LMXB) are radio emitting sources. Some (~20) of 
these radio emitting sources show a relativistic jet 
and are called microquasars (see Refs. 3 and 4 for 
reviews on microquasars). 

Some of these sources were observed with the 
first generation of TeV instruments. The source most 
widely observed was the HMXB Cygnus X-3 and al- 
though there were claims of its detection, it has not 
been considered as a positive one by the astronomical 
community because it has never been confirmed by 
the new generation of Cherenkov telescopes and the 
instrumentation at this epoch was limited. Another 
monitored source was the LMXB microquasar GRS 
1915+105, which was the first superluminal source 
detected in the Galaxy [5] . The HEGRA experiment 
detected a flux of the order of 0.25 Crab from GRS 



1915+105 during the period May- July 1996 when the 
source was in an active state [6] . This source has also 
been observed with Whipple, obtaining a 3.1c signifi- 
cance [7] . Although these sources have been observed 
with the new generation of Cherenkov telescopes and 
they remain undetected up to now, other binaries 
have been detected at TeV energies. 

Nowadays, the VHE 7-ray sky map contains four 
XBs, namely PSR B1259-63 [8], LS I +61 303 [9], 
LS 5039 [10] and Cygnus X-l [11]. Whereas the op- 
tical companion in all these sources is well known, a 
massive O or B star, the nature of the compact com- 
panion is unknown in two of them. While in the case 
of PSR B1259— 63 it is clear that the compact object 
is a rapidly spinning non-accreting NS and the TeV 
emission is powered by rotational energy [8, 12], and 
in the case of Cygnus X-l it is an accreting stellar- 
mass BH, in LS 5039 and LS I +61 303 the uncer- 
tainties in the determination of the inclination of the 
system [13, 14] precludes to fix the mass of the com- 
pact object and, therefore, to know if it is a BH or a 
NS. 

In this work, I report in Section 2 the main char- 
acteristics describing each of the four X-ray binaries 
detected at TeV energies up to now. In Section 3 
I present some scenarios for high energy (HE) and 
VHE 7-ray production. Finally, in Section 4 a brief 
discussion and summary is given. 

2. Binary TeV sources - BTV 

Here follow some of the main characteristics of these 
binary TeV sources (BTV). Some of these charac- 
teristics are shared while other are unique and make 
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each of these sources particular. Some properties are 
summarized in Tabic 1. 

2.1. PSR B 1259- 63 

PSR B1259-63 / SS 2883 is a binary system contain- 
ing a B2Vc donor and a 47.7 ms radio pulsar orbit- 
ing it every 3.4 years, in a very eccentric orbit with 
e=0.87. No radio pulses are observed when the NS is 
behind the circumstellar decretion disk (free-free ab- 
sorption). VHE gamma-rays are detected when the 
NS is close to periastron or crosses the disk [8] . The 
TeV light-curve shows significant variability and a 
behavior not predicted by previously available mod- 
els. The observed time-averaged energy spectrum can 
be fitted with a power law with a photon index IVhe 
= 2.7 ± 0.2 s tat i 0.2 sys 

Different models have been proposed recently to 
try to explain these observations. It is shown by 
Ref. 15 that the observed TeV light-curve can be ex- 
plained by the inverse Compton model under certain 
assumptions, although the large uncertainties, as well 
as the relatively narrow energy band of the available 
TeV data do not allow to put robust constraints on 
model parameters. In Ref. 17, a hadronic scenario is 
proposed where the TeV light-curve, and radio/X- 
ray light-curves, can be produced by the collisions of 
high energy protons accelerated by the pulsar with 
the circumstellar disk. In this scenario the VHE 7- 
rays are produced in the decays of secondary ir°, 
while radio and X-ray emission are synchrotron and 
inverse Compton emission produced by low-energy 
(< 100 MeV) electrons from the decays of secondary 
7r ± . Coordinated observations in the keV, MeV, GeV 
and TeV bands are very important to constrain the 
models, especially at the epochs close to the perias- 
tron and disk passages. 

2.2. LS I + 61 303 

LS I +61 303 is a HMXB system that shows peri- 
odic non-thermal radio outbursts on average every 
P or b=26.4960 d [18]. The system is composed of a 
rapidly rotating early type BOVe star with a stable 
equatorial decretion disk and mass loss, and a com- 
pact object with a mass between 1 and 4 M Q orbit- 
ing it every ^26.5 d in an eccentric orbit [13]. The 
periastron takes place at phase 0.23±0.02. Spectral 
line radio observations give a distance of 2.0±0.2 kpc 
[19]. The maximum of the radio outbursts varies be- 



tween phase 0.45 and 0.95. X-ray outbursts, start- 
ing around phase 0.4 and lasting up to phase 0.6, 
have also been detected (see Ref. 20 and references 
therein). Orbital X-ray periodicity has also been 
found using RXTE/ASM data [21], which currently 
reveal a broad maximum covering phases 0.4-0.6. 
Similar results have recently been obtained at higher 
energies with INTEGRAL data [22]. The X-ray spec- 
trum can be fitted with an absorbed power-law with 
values in the ranges iVn=0.45-0.65x 10 22 cm' 2 and 
r=1.5-1.9 [23, 24]. Recent Chandra observations 
have reported possible evidence of an extended X- 
ray emission from the source [25] . 

LS I +61 303 is also spatially coincident with 
3EG J0241+6103, a HE (> 100 MeV) source de- 
tected by EGRET [26]. The MAGIC measure- 
ments showed that the VHE 7-ray emission from 
LS I +61 303 is variable and possibly periodic, with 
maximum flux ~16% that of the Crab Nebula de- 
tected around phase 0.6 [9] . The VHE spectrum is fit- 
ted reasonably well by a power law: dN/ (dAdtdE) = 

(2.7 ± 0.4 ± 0.8) X l -12£(-2.6±0.2±0.2) cm -2 g -l 

TeV -1 . A preliminary analysis of new data, which 
confirms the former results, has been presented by 
Ref. 27, and is shown in Fig. 1. In addition, the or- 
bital TeV variability has been confirmed by the VER- 
ITAS Collaboration [28]. 

The discovery of apparently relativistic radio 
emitting jets in LS I +61 303 allowed to classify 
this XB as a microquasar [29]. This extended jet- like 
radio structure, of angular extension of 10-50 mil- 
liarcseconds, seemed to be precessing. New and re- 
cent high resolution radio images obtained with the 
VLBA during a full orbital cycle show an elongated 
structure changing its orientation at different orbital 
phases [30]. This behaviour accommodates well in 
a model based in the interaction of a pulsar wind 
with the dense equatorial wind from the Be star [31]. 
However, according to Ref. 32, this scenario presents 
some difficulties to explain the GeV luminosity of 
this source because the required spin-down power of 
the putative pulsar should be so high that the wind 
of the Be star could not collimate the radio emitting 
particles. 

2.3. LS 5039 

LS 5039 is a HMXB system composed of a compact 
object of unknown mass, 1.4-5 M Q depending on 
the binary system inclination, that orbits an 06.5 
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Table 1. The four X-ray binaries that are TeV emitters. 



Parameters 


PSR B1259-63 


LS I +61 303 


LS 5039 


Cygnus X-l 


System Type 


B2Ve+NS 


B0Ve+NS? 


06.5V±BH? 


09.7Iab+ BH 


Distance (kpc) 


1.5 


2.0±0.2 


2.5 


2.2±0.2 


Orbital Period (d) 


1237 


26.5 


3.9 


5.6 


-^compact (Mq) 


1.4 


1-4 


1.4-5 


20±5 


Eccentricity 


0.87 


0.72 


0.35 


~ 


Inclination 


36 


30 ±20 


20? 


33 ± 5 


Periastron-apastron (AU) 


0.7-10 


0.1-0.7 


0.1-0.2 


0.2 


Activity Radio 


Periodic (48 ms and 3.4 yr) 


Periodic (26.5 d and 4 yr) 


persistent 


persistent 


Radio Structure (AU) 


£ 2000 


Jet-like (10-700) 


Jet (10-10 3 ) 


Jet (40) ± Ring 


iradio(O.l-lOOGHz) (ergs -1 ) 


(0.02 - 0.3) x 10 31 M 


(1 - 17) x 10 31 


1 x 10 31 


0.3 X 10 31 


^X(l-lOkeV) ( er S s_1 ) 


(0.3 - 6) x 10 33 


(3 - 9) x 10 33 


(5 - 50) x 10 33 


1 x 10 37 


Lvhe (erg s _1 ) 


2.3 x 10 33 M 


8 x 10 33 < a ) 


7.8 x 10 33( » 


12 x 10 33 < a ) 


TvHE 


2.7 ±0.2 


2.6 ±0.2 


2.06 ±0.05 


3.2 ±0.6 


Remarks 




3EG J0241+6103 


3EG J1824-1514 




Unpulsed radio emission; 


(a) 0.2 <E < 10 TeV; (b) 


Time averaged luminosity. 
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Fig. 1. MAGIC measured fluxes from LS I +61 303 as a 
function of orbital phase [27]. 

main sequence star every 3.9 d and is located at 
2.5 kpc [14]. Several 5 GHz VLBI observations of 
the source have been conducted in the past with the 
VLBA+VLA, EVN and MERLIN, always showing a 
central core and bipolar radio jets [33, 34]. This led 
the authors to suggest that the system was a new 
microquasar in the Galaxy. LS 5039 was proposed to 
be the counterpart of an unidentified EGRET source 
at HE 7-rays [33] . This has been confirmed after the 
detection in VHE gamma-rays (TeV energies) by the 
HESS Cherenkov telescope [10], showing also a 3.9 d 
periodicity (see Fig. 2) [35]. 

The differential photon energy spectrum is vari- 



able with orbital phase. During the phases of the 
compact object inferior conjunction the spectrum 
is consistent with a hard power-law where Tvhe = 
1.85 ± 0.06 s tat ± 0.1 svs t with exponential cutoff at Eq 
= 8.7 ±2.0 TeV. At the superior conjunction phases, 
the spectrum is consistent with a relatively steep 
(rvHE = 2.53 ± 0.07 sta t ± 0.1 sys t) pure power-law 
(0.2 to 10 TeV). The HE/ VHE emission is basically 
interpreted as the result of inverse Compton upscat- 
tering of stellar UV photons by relativistic electrons. 
However, the nature of the acceleration mechanism 
that powers the HE/VHE radiation is unclear, and 
two excluding scenarios have been proposed. In the 
first one electrons are accelerated in the jet of a mi- 
croquasar powered by accretion [36]. In the second 
one they are accelerated in the shock between the 
relativistic wind of a non-accreting pulsar and the 
wind of the stellar companion [31]. In the pulsar 
scenario, a new model assuming two different sets 
of power-law spectral parameters of the interacting 
primary leptons has been developed to explain the 
HESS phenomenology [37]. Some important clues to 
understand the physics behind the TeV emission in 
LS 5039 have been recently presented in a scenario- 
independent work [16]. 

VLBI observations tracing the motion of the ra- 
dio core and characterizing the changing morphology 
of LS 5039 during a complete orbital cycle would al- 
low to unveil the nature of the particle acceleration 



March 7, 2008 14:18 



WSPC - Proceedings Trim Size: llin x 8.5in paredes 



1 



r 5 




Fig. 2. HESS integral 7-ray flux of LS 5039 as a function of 
orbital phase [35]. 



mechanism and constrain some physical parameters, 
like the relativistic jet bulk velocity and inclination 
in the microquasar scenario or the binary system 
inclination in the non-accreting pulsar scenario. In 
a recent paper, [38] new VLBA radio observations 
of LS 5039 are presented. The imaging results (see 
Fig. 3) show a changing morphology between two 
images obtained five days apart, whereas the astro- 
metric results are not conclusive. This precludes to 
exclude clearly any of the two scenarios. While the 
observed changes in morphology cannot be explained 
easily with a simple and shock-less microquasar sce- 
nario, the young non-accreting pulsar scenario re- 
quires the inclination of the binary system to be very 
close to the upper limit imposed by the absence of 
X-ray eclipses. 

2.4. Cygnus X-l 

Cygnus X-l is the first binary system where dynamic 
evidence for a BH was found [39]. According to the 
most recent estimates, the BH mass is 20±5 Mq, 
while the 09.7 lab supergiant companion has a mass 
of 40±10 M Q [40]. The orbit of the system is cir- 
cular, with a period of 5.6 days and an inclination 
of 33±5° [39]. Located at 2.15±0.20 kpc (3cr error; 
see Ref. 40 and references therein), Cygnus X-I is 
the brightest persistent HMXB in the Galaxy, radi- 
ating a maximum X-ray luminosity of a few times 
10 37 erg s _1 in the 1-10 keV range. 

The source displays the typical low/hard and 
high/soft states of accreting BH binaries, spend- 
ing most of the time, currently about 65%, in the 
low/hard X-ray state [41]. Steady compact jets are 



produced in BH binaries in this state, when the in- 
ner radius of the disk is thought to be truncated, 
while in the high/soft state the jet is quenched [42]. 
This is also the case for Cygnus X-l, which displays 
a ~15 mjy and flat spectrum relativistic compact 
(and one-sided) jet (v > 0.6c) during the low/hard 
state [43], transient relativistic jets (v ^ 0.3c) dur- 
ing state transitions [44], while no radio emission is 
detected during the high/soft state. 

Arc-minute extended radio emission around 
Cygnus X-l was found [45] using the VLA. Their 
disposition reminded of an elliptical ring-like shell 
with Cygnus X-l offset from the center. Later, as 
reported in Ref. 46, such structure was recognised 
as a jet-blown ring around Cygnus X-l. In analogy 
with extragalactic jet sources, the ring could be the 
result of a strong shock that develops at the loca- 
tion where the pressure exerted by the collimated 
jet, detected at milliarcsec scales, is balanced by the 
ISM. If this pressure were solely due to a magnetized 
plasma in equipartition then we should expect a lobe 
synchrotron surface brightness much higher than the 
observed upper limit. This suggest that most of the 
energy is stored in non-radiating particles, presum- 
ably baryons. 

The instrument COMPTEL, on board the 
CGRO, did detect Cygnus X-l in the 1-30 MeV 
range several times [47] . Unfortunately, EGRET per- 
formed rather few observations of Cygnus X-l and 
did no detect the source. Only a quite loose upper 
limit to the ^100 MeV flux from Cygnus X-l is avail- 
able [48]. 

MAGIC observed Cygnus X-l in 2006 obtaining 
evidence of 7-ray signal with a significance of 4.9 
a (4.1 a after trial correction) [11]. The signal was 
variable and extending in a short time interval (~ 
80 minutes), as can be seen in Fig. 4. The measured 
excess is compatible with a point-like source at the 
position of Cygnus X-l and excludes the nearby radio 
nebula powered by its relativistic jet. 

The differential energy spectrum is well fitted by 
a power law given as dN / {dAdtdE) — (2.3 ± 0.6) x 
10~ 12 (£;/1 TcV)- 3 - 2±0 ' 6 . The measured excess was 
observed at phase 0.91, very near the inferior con- 
junction of the optical companion and superior con- 
junction of the compact object [49]. According to ex- 
isting models dealing with photon-photon absorption 
and cascading, there should not be detectable TeV 
emission at these orbital phases if its origin is close 
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Fig. 3. VLBA maps of LS 5039 at different orbital phases [38]. 
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Fig. 4. MAGIC measured fluxes from Cygnus X-l as a func- 
tion of time [11], The minor thicks in the X-axis are separated 
five days. 

to the compact object [50] (behind the companion 
star in this configuration). 

Interestingly, this detection occurred at the time 
when a flare was detected [51] at hard X-rays by 
INTEGRAL (at a level of about 1.5 Crab (20- 
40 keV) and 1.8 Crab (40-80 keV)) and Swift/BAT. 
The MAGIC detection occurred during a particularly 
bright and hard X-ray flare that took place within a 
prolonged low/hard state of Cygnus X-l. The TeV 
peak appears to precede a hard X-ray peak, while 
there is no particular change in soft X-rays. 

3. Scenarios for HE/VHE 7-ray 
production 

The observational results in the field of point-like 
TeV emission from galactic sources discussed above 
show that very efficient particle acceleration takes 
place in these compact sources. Photons have been 
observed up to energies as high as 30 TeV. 

Different works published during the last years 
have been based in the microquasar scenario [36, 52, 



53]. In these models, it has been proposed that the 
HE and VHE 7-ray emission of XBs with a jet (mi- 
croquasars) is due to inverse Compton up-scattering 
of stellar UV photons from the companion star by 
the relativistic electrons of the jets. 

Alternative models based on hadronic processes 
within the jets of a microquasar, or in the interac- 
tion between the relativistic jets and the interstellar 
medium (ISM), have also been proposed [54]. Finally, 
a different scenario for the origin of HE and VHE 7- 
rays in these sources is the region of interaction be- 
tween the relativistic wind of a young non-accreting 
neutron star and the wind of the stellar compan- 
ion [31]. 

These models predict different multi-wavelength 
behaviors that can be probed observationally. In par- 
ticular, very-high resolution interferometric radio ob- 
servations (VLBI) can provide a final answer to this 
open question. 

4. Discussion and Conclusions 

Up to now, 4 HMXBs but no LMXBs have been de- 
tected at TeV energies. This fact points the impor- 
tance of having a bright companion (O or B star) as 
source of seed photons for the IC emission and target 
nuclei for hadronic interactions. The compact object 
in these binaries is in one case a BH and in another 
case a non-accreting NS, whereas there are two cases 
where the nature of the compact object is unknown. 

Periodic TeV emission is present in one system, 
LS 5039, and likely in other two, LS I +61 303 and 
PSR B1259-63. In Cygnus X-l the detected TeV 
emission comes from a flaring episode and the de- 
tection of a possible steady VHE flux is below the 
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present IACT's sensitivity. The TeV flaring activ- 
ity observed in Cygnus X-l might not be a unique 
and rare case. When looking at the integral 7-ray 
flux of LS 5039 (see Fig. 2) it seems that at around 
phase 0.8 there is some flare activity superposed 
to the periodic-regular light-curve. In the case of 
LS I +61 303 in Fig. 1 we can see the expected max- 
imum at phase 0.6 and, in addition, a flaring activity 
at phase 0.8. This phenomenology seems to suggest 
that flare TeV emission might be common in XBs, al- 
though more flare detections are needed to establish 
this kind of TeV activity. 

The increase of observational capability at VHE 
produced by VERITAS and CANGAROO III and 
the gain in sensitivity of HESS II and MAGIC II, 
will allow to detect new XBs and establish if the flar- 
ing activity is a common phenomena in these type 
of sources. All of this will bring more constraints 
to the physics of these systems. In addition, multi- 
wavelength (multi-particle) campaigns will also be 
very important to gain knowledge about these new 
TeV sources. 
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